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In this paper we discuss a technique for calculating moments of polydisperse materials in terms of 
concentration readings along the cell. The proposed method minimizes dependence on data from the 
end points where they may be unreliable. An analysis is given of the errors involved in the use of the 
proposed method when the underlying molecular weight distribution is the Schulz distribution or the 
lognormal. 
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One of the primary functions of a sedimentation 
equilibrium experiment is to measure the molecular 
weight of the solute, and in the case of a polydisperse 
solute, obtain information about the molecular weight 
distribution. In the latter case, the information is in 
the form of the first several moments of the distribu- 
tion. The commonly used methods of data analysis 
derive these moments from the values of the concentra- 
tion and its spatial derivatives evaluated at the end 
points of the solution column [l]. 3 Methods of eval- 
uating the moments using a point nearer the center of 
the solution column have been described by Fujita [2] 
and Adams [3j. However, these techniques require 
data over the entire range from zero to infinite centrif- 
ugal field. Thus far, no practical test of these latter 
methods has appeared in the literature. 

Recently it has been suggested [4] that improved 
accuracy could be obtained for the moments and hence 
the molecular weight averages if measurements were 
made as a function of centrifugal field as the field 
approaches zero. This latter method again requires 
some data to be obtained by extrapolation to the menis- 
cus and other data to be obtained by extrapolation to 
the cell bottom or from a point near the center of the 
cell. The purpose of this paper is to present a general 
method of analyzing experimental data in which the 
end points play a less important role and advantage is 
taken of the more accurate data obtainable elsewhere 
in the cell. It will also be shown that the treatment of 
Osterhoudt and Williams [4] represents a special case 
of the general treatment presented below. In addition 
to presenting the method we shall analyze possible 
errors in the method when the underlying molecular 
weight distribution is the Schulz [5] distribution or the 
lognormal distribution [1]. 
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Method of analysis. In the case of an ideal, non- 
compressible solution the radial concentration distribu- 
tion of a single solute species d(Q of molecular weight 
M i , is given by [2] 



c? "" l-exp(-\M,) ' 
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where cfis the original concentration before sedimenta- 
tion and where £ is the reduced radial variable given by 
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The distance of the meniscus and bottom of the solu- 
tion column from the center of rotation are given by 
a and b respectively, and r is any arbitrary intermediate 
position. The quantity A, is defined by 
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where the other symbols have their usual definitions [2]. 
For a solute which is not monodisperse but has a 
continuous distribution of molecular weights given by 
some function /(M), the radial distribution of the con- 
centration at sedimentation equilibrium is given by 
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where 0(Q is defined by this equation. 

One observes that the form of the function given by 
the right-hand side of eq (1) and contained in the 
integrand of eq (4) is the same as that of the generating 
function for the Bernoulli polynomials [6], provided 
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that KM ^ 277. Thus eq (4) can be rewritten 



where a can be chosen arbitrarily. By choosing a— Va 
we find 
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where i>J is the truncated moment 
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or choosing a = 0, 
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M$, the truncated weight average molecular weight 



Values of the moments are derived from these relations 
by taking measurements at several values of X and 



and so on. Some properties of the polynomials B„(Q extrapolating to the dependence at X = 0. 



are given in reference 6. 

The basis for the method suggested by Osterhoudt 
and Williams is the replacement of 0(£) by 0*(£), 
defined by 
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i.e., the integral in eq (5) is assumed to be negligible. 
With this definition, the identity d B„(Q I d£.= n B„-x(Q, 
and the particular values B (0) = 1, B,(0) = — 1/2, 
B-i0) = l/6, B :i (0) = 0. B 4 (0)=-l/30, B 5 (0)=0, ... 
we readily obtain the result of Osterhoudt and Williams 
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The remaining formulas in reference 4 for £ = 1/2 and 
1 are derived in similar fashion. However, we can 
also derive other identities that allow us to use any 
£ values between zero and one thus permitting the use 
of points that obviate extrapolation and are therefore 
more reliable. For example, we have 
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It is possible to extend these considerations so that 
a series is obtained, the first term of which is propor- 
tional to \ r v r , and in which any number of terms in 
v*, v%, . . . , v*_ x , v? +1 , . . . Vr + m have zero coeffi- 
cients, by properly choosing points £i, £2, • . • , (,r + m, 
and coefficients a u a 2 , . . . a r + m, in linear com- 
bination: 



/ r *(£l,£ 2 , • • • ,£r+m) = C*10*(£l) 

+ « 2 0*(£ 2 ) + . 

Given the set of points £1, £2, 
chosen so that 
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One of the constants, a u can be set equal to 1 and the 
ratios otj/au are determined from eq (13). The points 
£/ can be chosen arbitrarily except for the condition 
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in order that the coefficient of \ r v * be nonzero. As an 
example we can calculate a formula in which the first 
nonzero is proportional to k 3 v* , and the second XV*, 
by choosing 

£i = l/5, £ 2 = 2/5, £.3 = 3/5, £4 = 4/5. 

Using the procedure outlined above we find that 
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Error analysis. So far we have made the tacit 
assumption that the observed value 0(£) is identical to 
0*(£) so that v* can be identified with the desired value 
v n . It clearly is not, so that some error analysis is 
required to set bounds on the validity of the method. 
Two types of error require investigation. The first 
is the error in using the observed 0(£) rather than the 
required 0*(£), and the second is in the calculation of 
values of v% rather than v n . 

Let us denote by e(£) the difference 
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The absoulte error incurred by using values of 0*(£) 
rather than 0(£) can be bounded as follows: 
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so that the difference between the observed F and the 
desired F* is proportional to k(v\ — v*). Theoretically 
one can decrease this bound by choosing the £j close 
to 1 but this method is limited by observational errors 
the meniscus. It is interesting to note that if the 
moments are to be calculated in terms of parameters 
relating to the derivative 0'(£), as is the case in Oster- 
houdt and Williams' paper, the bound corresponding 
to eq (17) contains \(yi — v*) in place of k(v x —v*) and 
is therefore larger. 

For the purpose of illustrating errors in terms of 
physical quantities we will derive explicitly the for- 
mulas for \(i>i — P*) and Jn=l—(v%lvn) for two com- 
mon polymer distributions; the Schulz [5] and the 
lognormal [1J. We consider first the Schulz distribution 
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where the two adjustable parameters a and p can be 
expressed in terms of M w and M z , the weight and z 
average molecular weights as 

a = (M z -M w )~ l 



We therefore see that the difference depends on the 
parameters p and k/a where 



k/a = k(M z -M w ). 
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Experiments of the type discussed here can be ar- 
ranged so that kM w ~ 1. Furthermore, the Schulz 
distribution is a sensible one for polymers only for 
p positive or 2 > MJM W hence kM z ^ 2 and XI a is 
at most equal to 2. For pk/(27ra) a small number we 
can approximate the integral in eq (20) by 
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In figure 1 we have plotted some representative 
values of — logi X(^i — vf) as a function of kM w . As 
can be seen from the figure 6(Q and 0*(£) are experi- 
mentally indistinguishable when kM u is less than 1 
except when M Z \M W is greater than 1.75. Even in 
that case if \M W can be set less than 0.5, 0(£) and 0*(£) 
are experimentally indistinguishable. 



To estimate the accuracy with which 
mates to v n we calculate the ratios 
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In figure 2 we have plotted — lognJi and — logi(x/ 2 
as a function of kM w for M g /M w =1.75. It is to be 
noted that both v\ and v 2 can theoretically be deter- 
mined to within an error of about 1% with the present 
method, provided that kM w is less than 1. 

Another distribution useful in polymer work is the 
lognormal [1] 



f(M)- 



1 



o-MV2tt 



exp 



1 
'2a 2 



M_ 
Mo. 



(25) 



_ 2M W -M Z 
P M z -M w 

For this distribution X.(^i — v* ) can be written 
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The two adjustable parameters M and cr, are related 
to M w and M z by 
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straightforward calculation suffices to show that 
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k n {v n — v*) can be expressed as 
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For odd n these can be written in terms of complemen- 
tary error functions and exponentials, while for even 
n they can be written in terms of exponentials. The 
J n defined in eq (24) are expressible in the form 

Jn= | u n ~ x e- u ^du\ \ u n ~^e- u ^du. (28) 
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Specifically, the first three 7's are 
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FIGURE 1. Graphs of —log 10 X(^i — ^*) as a function of XM W for 
different values ofMJM w for a Schulz distribution. 
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FIGURE 2. Graphs of —Log 10 Ji and —log 10 J 2 as a function of \M„ 
for Schulz distribution. 
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where erfc(. 
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These formulas can be simplified by noting that 
practical values of \M (T are of the order of 1 or less. 
As we have already noted, KM W can be made of the 
order of 1, and M W \M Z < 1 so that XM < kM w . Fur- 
ther, for the cases where the ratios (M z /M w ) are less 

2tt 
than approximately 10 and cr < 1.6, the ratio -rjz — is 
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greater than 3 and we can approximate the complemen- 
tary error function by 

r ( 2tt \ KMo<r ( 2tt 2 \ 

erfc \xm^) ~ -*r exp Vm^V (30) 

which can be inserted in the expressions for A(i>i — vf) 
and Ji, J 2, J3. This procedure yields 
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FIGURE 3. Graphs of —log, X(^i~ v*) as a function of MJM W 
for lognormal distribution. 



Figure 3 shows curves of — logioM^i — rf) as a function 
of MJM W for \M, r = 3 and 4. It can be seen that 
kv x is very close to \v^ for kM w = 3 and all MJM W . 



Since \M W can be kept in the neighborhood of 1 the 
present technique of deriving moments from the con- 
centration can be justified theoretically for all M Z \M W 
provided that the lognormal distribution is a suitable 
representation of the weight distribution. Calculation 
of the 7's serves only to confirm this observation. It 
is of some interest to note that in contrast to the results 
for the Schulz distribution, the } n decrease when 
M z jM w increases. 

The two distributions assumed here are for illustra- 
tive purposes only. It is probably safe to say that if 
\M W ^ 1 the errors made in the mathematical assump- 
tions are negligible compared to the experimental 
errors. 
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